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Q>^ I Isolated electrons resting above a helium surface are predicted to have a bound 
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' spectrum corresponding to a one-dimensional hydrogen atom. But in fact, the 
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observed spectrum is closer to that of a quantum-defect atom. Such a model is 
discussed and solved in analytic closed form. 
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Some time ago the prediction was made that an isolated electron resting on a helium (or 
some certain other) surface should have a bound-state spectrum in the vertical direction 
The idea is that the electron induces an image charge in the helium, producing a potential on 
the electron of 

= +00, X < 0, (2) 

where e is the dielectric constant [^. For helium it is |5|, ^ 

e = 1.05723, Z = 0.0069547. (3) 

The spectrum should thus be similar to that of a (weakly-coupled) one-dimensional hydrogen 
atom. This phenomena has been observed (See for a current review and |10| for a 

proposed application to quantum computing.) 



Consider the one-dimensional Schrodinger equation of this system |11]: 



y 2m dx^ 
Making the changes of variables 



1pn{x) = En1pn{x). (4) 



one obtains 

Observe that the helium-surface "Rydberg" and "Bohr radius" have values 

So = Z'^Roo = 0.658086 meV = 159.123 GHz, (8) 
bo = 2x0 = ao/Z = 76.01 A. (9) 

By techniques similar to those used to obtain the solutions for the 3-dimensional hydrogen 
atom, one can obtain the normalized eigensolutions |12, ^] 



Mz) = pn^xo]"'/' zn exp[-z„/2] L^^l,{zn). (10) 



This agrees with particular n = 1, 2, 3 wave functions in the hterature [p^ ]. 

In Eq. (^) we have used the generalized Laguerre polynomials commonly found in the 
modern mathematical physics literature [p^]: 



A:=0 \ 



k\ 



(11) 



In Eq. (pA|), the generalized binomial symbol ( ^ ) means r(a + l)/[r(a — 6+ l)r(6+ 1)]. Also, 
Ln\x) = Ln{x), where Ln{x) are the ordinary Laguerre polynomials normalized to unity at 
zero: Ln{0) = 1. These polynomials were used instead of the associated Laguerre polynomials 



often defined, for Coulomb wave functions [|16[, as 



Liix) 



d^Lnix) d^ 



dx^ 



dxi 



dx^ 



(12) 



Here, Ln{x) = (n!)L„(x), the ordinary Laguerre polynomials normalized to L„(0) = (n!). Eq. 



(12) can be confusing, since this definition only holds for integer j. Contrariwise, Eq. ( [Tl| ) is 
defined for arbitrary a. (This will be very important in the following.) When a = j, an integer, 
the connection between the two forms is 



Li^^ix) = {-iy[{n + jy.] L^\x). 



(13) 



The experiments obtain transition energies from excited states to the ground state: 



1 



An — \Ei \ — \En\ — £"0 1 :t . 



(14) 



However, the experiments do not yield exact Balmer energies. The A„ are all of order 7 GHz 
too large [^-[0]. This is like a quantum defect, since if 

£0 



then 



En En* 



(n*)2' 



n* = n — (5, 



(15) 



An- 



l-E'l* I — \En* I — £q 



(1-5)2 {n-5f 



£0 



1- J-) +25(1-^) + 36^(1-^]+.. . 



(16) 
(17) 



The correction term, £q 2 5 [1 — 1/n^], varies by only ~ 10 % as n varies from 2 to oo [17|. 
Just fitting the 2* 1* and 3* 1* transition energies to this formula yields, Sq = 158.4 
GHz and 6 = 0.0237 or an increase in A of about 7.8 GHz. In other words, this is like a 
one-dimensional Rydberg atom. 



Elsewhere |18||-[pi|, inspired by supersymmetry [22|, it was shown how one can obtain exact, 
analytic, one-particle wave functions for real Rydberg atoms yielding the correct eigenenergies. 
This also yielded: (a) transition matrix elements in agreement with experiment and complicated 



many-body calculations ||l^; (b) good fine-structure splittings |19|; and (c) Stark splittings 



whose crossing/anti-crossing patterns agree with experiment [20|. 

The mathematical key to this success is the fact that for proper solutions of the (radial) 
hydrogen- atom equation one does not really need that / = (integer) and n = (integer). One only 
needs that (n — /) =(integer), the two separately not having to be integers. That is, the factor 
/(/ + 1) in the effective 1/r^ potential term need not have I be an integer for a finite-order 
polynomial radial solution to exist. This is where the L^\x) become of use p^]. 

Applying this idea to the present case, we phenomenologically propose for x > that V{x) 
becomes 



V{x) 



Ze^ , (-5)[{-6) + l] 



+ 



X 



2m 



X > 0. 



x^ 



(18) 



Then the exact eigenenergies are given by Eq. (^) and the exact wave functions are 



Nn* 4* exp[-z„V2] L'^^_t'\zn'), 
1 r(n) 



(19) 
(20) 



.2(n*)2xo r(n+ 1 - 2(5)J n'^xo 

Thus, we have an exact analytic solution to the problem. We can also analytically calculate the 
expectation values [^] as double sums of Gamma functions |^. In particular Jl3[ , 

{x)n* = xq hn^ - 5(6n - 1 - 25)] . (21) 



When 5 = this reduces to the standard result. Also 



(l*|x|n*) 



r(n + l-25)r(n) 
r(2 - 26) 



.1/2 



n-1 

E 

k=0 



{-gf {k + 2,-25){k + 2-26) 



k\ 



{n-l-k) 
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2 ■ 1* 

n* + 1* 



(22) 



Setting n* to 1* in Eqs. ( pi| ) and (22), makes then equal. 

Unfortunately, this model does not resolve the physical problem of how one realistically cuts 
off the unphysical, negatively-infinite potential at the origin ||2^. In fact, this solution makes 
the problem slightly worse: at the origin the potential now goes to negative infinity as —l/x^. 
If the experimental quantum defect had been of opposite sign, then the added potential would 
have been positive, like an angular momentum barrier, making the states less bound. This also 
would have " realistically" modeled the positive work function at the surface of about 1 eV ||8|. 



There are inverse methods for generating inequivalent isospectral Hamiltonians |27]-[30|. 
What, in principle, would be an isospectral Hamiltonian with the desired physical properties 
is one with an added potential that: (i) goes, at the origin, to plus infinity at least slightly 
faster than [5(1 — 5)\/z^ , (ii) becomes negative for larger z, and (iii) goes to zero at infinity 
from below. 

A first examination of the above inverse methods |27|- [pO| found potentials with the last two 
properties, but not the first. These potentials go to zero at the origin. An example is [27, IsO] 



V2{Z) 



(1 



*\2 



26 



l]Y + Y' 



X > 0, 



Y 



exp[-zi*] 4-2^ 



r(3 - 26, zi*)-R 

where R (which can be chosen to be —2) and r(a, z) (the incomplete F function) are 
7 + 1 7 + 1^ 



(23) 
(24) 



R = 



7 (1*) Nl 7 



-r(3-25). 



T{a,z) 



dy y 



(25) 



and 7 is a dependent constant useful below. Taking units of xq = 1, the orthonormal eigen- 
functions are 



Xn*{z) 



Kiz,y) 



tpn*{z)+ dy K{z,y) %IJn*{y), n > 1, 
Jo 

1 X exp[-zp/2] (z\7^) exp[-2/i./2] (y\7^ 



1* 



r(3 - 26, zi*)-R 



(26) 



(27) 



The exception, with normahzation 7/(7 + 1), is 

7I/2 )t{3-26,z,.)-R- ^^^^ 

There may well be analytic isospectral Hamiltonians with all the desired properties. But to 

determine their existence requires more investigation. 
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